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Recently an ensemble of nuclear spins in a quantum dot have been proposed as a long-lived 
quantum memory. A quantum state of an electron spin in the dot can be faithfully transfered into 
nuclear spins through controlled hyperfine coupling. Here we study the decoherence of this memory 
due to nuclear spin dipolar coupling and inhomogeneous hyperfme interaction during the storage 
period. We calculated the maximum fidelity of writing, storing and reading operations. Our results 
show that nuclear spin dynamics can severely limits the performance of the proposed device for 
quantum information processing and storage based on nuclear spins. 



INTRODUCTION 

Electronic and nuclear spins in semiconductor nanos- 
tructures have been studied intensively because of their 
potential applications in quantum computation and spin- 
tronics. In a series of proposals of solid state quantum 
computers, electron or/and nuclear spins have been sug- 
gested as qubits |l|,|2|, the information carriers. Hyper- 
fine coupling, the main interaction between an electron 
and the surrounding nuclei in semiconductors, can be uti- 
lized as a tool to polarize the nuclei by optical pumping 
P| or spin-polarized transport [j, la, la]- On the other 
hand, it could also be a major decoherence channel for 
electron spins yj B 

In a recent study, long-lived quantum memory of the 
electron spin in both quantum dots 11(1 and quantum 
Hall semiconductor nanostructures |llLll2j . with the help 
of hyperfine interaction, is proposed using highly polar- 
ized nuclei in these structures. The idea is to transfer the 
electron spin state to nuclear spins coherently by tuning 
to the resonance condition of the nuclei-electron hyper- 
fine flip-flop term; this process is referred to as informa- 
tion writing. Since nuclear spins have a long relaxation 
time T\ (up to minutes), it is hoped that the informa- 
tion of electron spin can be stored in nuclei for a much 
longer period of time. The reverse process (reading) can 
be activated on demand later by tuning back to the res- 
onance point so as to transfer the collective nuclear spin 
state back to the electron spin. Partially polarized nuclei 
which uniformly couple to the electron have been shown 
not to depreciate the performance of the quantum mem- 
ory dramatically [lOj . Internal dynamics of nuclear spins 
has so-far been neglected. 

However, nuclear spins are not static. Their evolu- 
tion during the storage time, when their flip-flop with 
the electron is turned off by either changing the external 
magnetic field or simply removing the electron from the 
dot, will induce decoherence in the quantum memory. In 
this paper we study the influence of nuclear spin dynam- 
ics on the information transfer and storage. Two types of 
spin interactions of nuclei could be in effect, namely the 
dipole-dipole coupling and the inhomogeneous hyperfine 
coupling with the electron (effective field seen by nuclei). 



During the writing and reading processes, the dipolar in- 
teraction can be ignored, since the operation time r is 
determined by the hyperfine coupling constant (see be- 
low), so that rHd/fi is much less than 1, where Hd is the 
dipolar Hamiltonian. However, it could be critical when 
the storage time approach the nuclear transverse relax- 
ation time T2, which is usually in the order of 10~ 3 to 
10~ 4 second. 

In this paper we assume that initially the nuclear spins 
are fully polarized in a pure state \I, ■ ■ ■ , /), where I rep- 
resent the magnitude of nuclear spins, and mi means 
complete polarization. The system will undergo non- 
trivial evolution when a spin down electron is put into 
the dot. Once spin transfer is achieved at t w , we turn off 
the resonance condition and let the nuclear spin system 
evolves for a long time r s . Finally we read the informa- 
tion by putting back a spin up electron and turning on 
the resonance. Our objective then is to study how the 
maximal reading fidelity vary as a function of the storage 
time t,. 



COHERENT SPIN DYNAMICS OF THE 
EFFECTIVE TWO-LEVEL SYSTEM 

In the writing and reading stage, the electron-nuclei 
spin system can be approximated with a simple Hamil- 
tonian 
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The first two terms in H are the Zeeman energies of the 
electron spin and nuclear spins. Hh describes the hyper- 
fine interaction in which a single electron interacts with 
an ensemble of N nuclear spins where N m 10 4 — 10 5 for a 
small quantum dot; Ai = ^4|-0(Ri)| 2 is the local strength 
of hyperfine coupling and ip(Ri) is the electronic enve- 
lope wavefunction at the zth nucleus. In the Hamilto- 
nian H we neglect the nuclear dipole-dipole interaction 
as we have previously explained. The combined system 



of nuclear and electron spins can be taken as an effective 
two-level system, where the electron spin could be cither 
up or down. The time-dependent wavefunction reads 
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where |* ) = | ^;/,---,/) and |tf fc ) = | fr, /,..., (/ - 
l)fc, • ■ • , I). Substituting Eq. |2Jinto Schodingcr equation, 
ihd$(t)/dt = HV(t), we find 
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where <5 = w e + AI/K and 7 = NIuij. Under the reso- 
nance condition, the electron-nuclei flip-flops are greatly 
enhanced since energies arc conserved during these pro- 
cesses. The resonance condition can be obtained by 
equating the total energies before and after the scatter- 
ing. We arrive at 
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In deriving Eq. 0] we have assumed that the hypcrfinc 
coupling constant is A/N. Since N sa 10 5 3> 1 and 
Wj w 10~ 3 cu e -C W e , Eq. 0]can be reduced to £ = w e + 
AI/h « 0, i.e., the resonance condition can be satisfied 
by adjusting the external magnetic field. 

An exact solution of Eq. |21 was found using Laplace 
transformation pj. However, the non- trivial complex in- 
tegral (the inverse Laplace transformation) in the solu- 
tion makes it less transparent. In the following we will 
apply a different approach. For later convenience we de- 
fine 
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a and y\. This approximation is appropriate because 
the decoherence induced by the inhomogencous hypcrfine 
coupling is significant only when the time scale is as long 
as N/A for fully polarized nuclei pfl. The spin transfer 
time in which we are interested is much shorter. The 
argument will become clear in the following calculations. 
We consider two initial conditions corresponding to the 
writing and reading processes respectively. For the initial 
state with a spin down ^(0) = | JJ-; I, ■ ■ ■ ,1); a(0) = 1 and 
Pi = ■ ■ ■ = Pn = and 3/1 (0) = 0. The solution of Eq. 
and Eq. His 
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where il = ./2I^. A 2 - is the Rabi flopping frequency of 
the effective two-level system. Substituting Eq. [HJback 
into Eq. [3] one find the solution of P k is 
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We see that the system oscillates coherently between 



state |*o) and state £ fc A k \^ k )/y/J2j A y If we kee P 2/2 
in Eq. the solution will be different from Eq. [5] How- 
ever the corrections of the amplitude and the frequency 
for a arc ~ 1/vN and A/N respectively. In other words, 
the approximation Eq. [3 we made should be OK as long 
as the time scale is much smaller than N/A. If we let 
t w = 7r/Q, |a(tm)| 2 = so that complete spin trans- 
ferring is achieved. The condition t w <C N/A is consis- 
tent with our approximation. For the initial state with a 
spin up electron, *'(0)) = Efc#fc(0)l*fc); «'(°) = ° and 
2/i(0) = V ■ A?'A;(0)/^- We find the solution in this case 
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Using the resonance condition (5 = and Eq. [S] Eq. [31 

can be transformed into 
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Without further approximation Eq. El are equivalent to 
Eq. |3 However, if we write j/2 as 
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which is exact if the hypcrfinc coupling is uniform. The 
group of differential equations are now closed with only 



Suppose nuclear spin dynamics is frozen during the stor- 
age time, yi(0) = J2j A j^'j(°)/ h = Y. 3 A jPj{ T w)l^ and 
I a (7r/f2)| 2 = 1. However, if we allow the nuclear spins 
to evolve during the storage time, P k (0) will be differ- 
ent from p k (t w ), which will make \a (7r/f2)| 2 less than 1. 
In the next section, we shall investigate how the deco- 
herence induced by the nuclear spin dynamics affect the 
coherent electron spin transfer in the read out process. 



NUCLEAR SPIN DYNAMICS 

During the period of storage we can either turn off 
the resonance by varying the external magnetic field or 
simply remove the electron from the dot. In the first case 
the electron-nuclei flip-flops are largely suppressed; nuclei 



only feel the effective magnetic field generated by the 
term AiIi Z S z which acts like an inhomogeneous field for 
nuclei at different lattice sites. The effective Hamiltonian 
of nuclear spins is 
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Here Hd [lj| describe the secular terms of the dipolc- 
dipole interaction; fry = h 2 rf(l — 3cos 2 0y)i?~- 3 /4, where 
Rij is the distance between the nuclei at site i and site j 
and Oij is the angle between Ry and the external mag- 
netic field. The second term in Hd is the flip-flop term 
which introduces the nuclear spin diffusion 14]. If the nu- 
cleus at the center of the dot is flipped down at t = t w , 
then the down state of the nuclei will propagate to the 
nuclei at larger radius through their mutual flip-flops. 
After a sufficiently long period of time, the nucleus with 
flipped spin state may be at the very boundary of the 
dot. If one turns on the electron- nuclei resonance at this 
time, the injected spin-up electron has minimal chance to 
interact with the nucleus with flipped spin, because the 
probability of the electron at the boundary is negligible. 
This leads effectively to decoherence of the nuclear spin 
memory. 

Since the truncated dipolar Hamiltonian Hd conserve 
the total spin of all the nuclei in the quantum dot, the 
general time-dependent spin wavefunction can be ex- 
pressed in terms of a linear combination of \ipk) 
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which is a constant. The reason that we do not drop 
the second term on the right hand side of Eq. ^J al- 
though it is much less than 77, is that we are interested 
in the large time behavior (3> 7r/£7) where non-uniform 
hyperfine coupling may have significant effect. 

Even for such a simple situation, no analytical solution 
of Eq. EH can be found. Therefore we have to rely on 




FIG. 1: \/3(x, 50, 30)| 2 as a function of spatial coordinate x. 
The various curves in the figure represent the temporal pro- 
file of |/3| 2 at different times starting from t — 0. The time 
separation between every two neighboring curves is 1 millisec- 
ond. It has been assumed that L r — 25 and L z — 10 in the 
simulations. The external magnetic field is along the x di- 
rection. The lattice size which we use in these calculations is 
101 x 101 x 61. The coordinate system has been established 
so that the center of the dot is located at (50,50,30). 




FIG. 2: |/3(50, 50, 2;) | 2 as a function of spatial coordinate z. 
All the parameters used in the calculations are the same as 
those in Fig. Q 



some numerical solutions. For this purpose we perform 
the following transformation for j3k(t) 
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Eq. 1131 then becomes correspondingly 
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This set of transformed equations shows that the factor 
r] only contributes an overall phase factor for every ftk, 
so it will not affect our calculation of the probabilities. 
We set the initial condition as /3fc(0) = i^U-/ V A 2 , i.e., 
the system starts right after the spin of the electron is 
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FIG. 3: Fidelity as a function of the storage time. The solid 
(dashed) represents the result of calculations when the exter- 
nal magnetic field is along z (x) axis where there is no electron 
in the dot (no hyperfine coupling) . The dot-dashed line shows 
the result when the hyperfine coupling is present. In all these 
calculations we assume a cubic lattice with size 101 x 101 x 61. 



transferee! to the nuclei (see. Eq. |SJ). We assume in our 
calculations that 
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where a = 5.65 A is the lattice constant for GaAs. 
The coordinates x,y and z are in units of a. L r and 
L z are two constants determined by the geometry of 
the quantum dot. I = 3/2 for all the nuclei in GaAs; 
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The time scale for dipolar interac- 



tion is a 3 /rfh « 10 milliseconds. We use both leapfrog 
and iterative Crank- Nicolson [15j method to solve Eq. ^j] 
We find that the leapfrog scheme is better in maintain- 
ing the unitarity of the coherent evolution for the nuclear 
spins, the total probability is preserved in the evolution. 
We assume a closed boundary so that spin can not dif- 
fuse out of the dot. For simplicity we also assume that a 
single species of nuclei resides on a simple cubic lattice. 
Let's now consider the case where the electron is re- 
moved from the dot so that the evolution of the nuclear 
spins arc only governed by the flip-flop term (nuclear spin 
diffusion). In Fig. ^ an d Fig. [2] we show the probability 
|/3| 2 as a function of spatial coordinates x and z. The 
slightly different behaviors observed in Fig. ^ and Fig. 
121 are due to the sharper confinement along the z direc- 
tion; we have also integrate Eq. 1 161 assuming the external 
field is in the z direction. Similar results as Fig. Q] and [21 
are found. After time t s , the coherent spin dynamics of 
electron can be activated by sending an electron with up 
spin into the dot. To calculate the fidelity of the reading 
process, we use the definition given in Ref. [lOj 
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F = Tr{p(t r )(S z - -Y}, 
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p(t) = \^(t))(^f(t)\, where the time-dependent wavefunc- 
tion takes the form of Eq. ^| It is then found that 
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where r r = 7r/0. The density operator is defined as 



The calculated fidelity is shown in Fig. Elby the solid and 
dashed lines for different external magnetic field orienta- 
tions. We see that the fidelity drops dramatically when 
t s reaches a few milliseconds. The difference of these two 
curves can be understood as follows. Nuclear spin dif- 
fusion is stronger along the external field direction. If 
this direction is the same as z direction, which has larger 
confinement (larger gradient) of the electronic wavefunc- 
tion, the diffusion process is even faster. This explains 
that fidelity represented by the solid line (z) is less that 
that of the dashed line (x) at equal time r s . 

It appears therefore desirable to suppress the decoher- 
ence induced by spin diffusion with an inhomogeneous 
hyperfine coupling |8j, |lj| . This can be realized by leav- 
ing the electron in the dot during time r s . This suppres- 
sion is due to energy non-conservation in the flip-flop pro- 
cesses. However our result (the dot-dashed line in Fig. I3J) 
shows that the coherence can only be maintained within 
a few microseconds, even worse than when the electron 
is removed. We can explain this result by identifying 
another dccohcrence channel, the inhomogeneous hyper- 
fine coupling itself. The nucleus at different lattice sites 
have different precession frequencies due to their interac- 
tion with the electron, so the phases of the wavefunction 
Pk{t) will not be uniform throughout the dot (see Eq. 
I15JI . This induces destructive interference of the function 
J2k Ak/3k(t). The time scale of this decoherence is deter- 
mined by the hyperfine coupling constant which is three 
order of magnitudes larger than the dipolar constant. 



CONCLUSION 

In conclusion, we have studied nuclear-spin dynam- 
ics induced decoherence of the long-lived quantum mem- 
ory proposed in Ref. [1CJ]. We find that storage time 
is limited to a few milliseconds (without hyperfine cou- 
pling) and a few microseconds (with hyperfine coupling 
AiIi Z S z ). The performance of nuclear spin memory is 
limited by the two decoherence channels: nuclear spin 
dipolar coupling and inhomogeneous hyperfine coupling. 
We only consider the fully polarized nuclei. Partially po- 
larized nuclei give rise to further decoherence. To make 
the device work as a truly long-lived quantum memory, 
one has to improve the device by combining it with other 
techniques such as rcfocusing. 



* Based on work presented at the 2004 IEEE NTC Quan- 
tum Device Technology Workshop. 
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